O Shell Effects in Bubble Nuclei, Atomic Clusters, and 

q ■ Inhomogeneous Neutron Matter * 

^ ■ AUREL BULGAC 1 AND PlOTR MAGIERSKI 2 

O 

department of Physics, University of Washington, Seattle, WA 98195-1560, USA 
^ i institute of Physics, Warsaw University of Technology, Warsaw, POLAND 



> . 

■ We analyze the character of the shell effects/Casimir energy in inhomo- 

geneous fermion systems. We estimate magnitude of the shell effects and 
discuss their dependence on a number of physical parameters (geometry, 
fermion density, temperature). 
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The problem of relative arrangement of impurities inside a many fermion 
system is determined solely by the quantum effects, as the volume, surface 
or curvature terms in the liquid drop expansion of the total energy are not 
affected when impurities are moved around (if Coulomb effects are irrele- 
vant). The properties of the quantum systems that contain impurities of 
^ ■ various shapes have been studied in the case of the Bose-Einstein conden- 

sate [1]. There, due to the fact that all the particles have the same single 
particle wave function, the behavior of the impurity is rather obvious and 
one can easily show that it will be expelled from the condensate [1]. In 
the case of fermion systems the many-body wave function has a complex 
character due to the Pauli principle and the answer is not obvious. 

The additional energy associated with the relative arrangement of im- 
purities can be termed shell correction energy. In fact, there is no well 
established terminology for the energy corrections we are considering here, 
even though the problem has been addressed before to some extent by other 
authors. In the case of finite systems, the energy difference between the true 
binding energy and the liquid drop energy of a given system is typically re- 
ferred to as shell correction energy. In field theory a somewhat similar 
energy appears, due to various fluctuation induced effects and it is referred 
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to as the Casimir energy [2]: 

/oo 
dee[g(e,l) - g (e)}, (1) 
-oo 

where go{e) is the density of states per unit volume for the fields in the 
absence of any objects, g(e, 1) is the density of states per unit volume in the 
presence of some inhomogeneities and 1 is an ensemble of geometrical param- 
eters describing these objects and their relative geometrical arrangement. 
Thus the Casimir energy can be thought of as a measure of the fluctuations 
induced in the energy spectrum in the presence of various "obstacles". 

There are a number of situations when such inhomogeneities can be 
formed. In the nuclear physics it was suggested a long time ago that very 
heavy nuclei will tend to develop a hole inside in order to minimize the 
Coulomb energy [3, 4, 5] and similar objects can be created in the case of 
atomic clusters [6]. 

The creation of voids is also predicted to happen in the nuclear mat- 
ter at subnuclear densities. Apparently, an agreement has been reached in 
literature concerning the existence of the following chain of phase transi- 
tions as the density increases: nuclei — ► rods — > plates — > tubes — > bubbles 
— > uniform nuclear matter [7]. In these cases, fermions reside in a rather 
unusual mean-field, which is much deeper for bound than for unbound nu- 
cleons. Since the amplitude of the wave function in the semiclassical limit 
is proportional to the inverse square root of the local momentum, the single 
particle wave functions for the unbound states will have a small amplitude 
over the deep well. Hence the deep well will act almost like a hard wall (in 
most situations). For the same reasons the halo nuclei [8] can be thought 
of as some kind of bubbles as well. 

Let us consider first the problem of positioning of a spherical bubble 
inside a finite spherical Fermi system. The total energy can be expressed in 
the form: 

E(N) = E LD {N) + E shell {N) = e v N + e s N 2 / 3 + e c N^ + E shell (N), (2) 

where Eld is the smooth liquid drop part of the total energy and E s h e u is the 
quantum shell correction contribution to the total energy. As we mentioned 
above, once the bubble is formed its displacement will not affect either the 
volume, surface or curvature terms in the liquid drop expansion. Hence, 
classically moving bubble off-center costs no energy, if Coulomb energy is 
left aside for the moment. The closer investigation of the shell correction 
energy shows that it depends strongly on the position of the bubble. The 
most pronounced shell effects are predicted for the spherically symmetric 
system (bubble in the center) . In this case the system is classically integrable 
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and thus its quantum-mechanical counterpart will have large gaps in the 
spectra resulting in a large amplitude of the shell correction energy. The 
rapid fluctuation of the shell energy as a function of the fermion density 
indicate that for some numbers of fermions it would be energetically more 
favourable to expell the bubble off-center. For a finite eccentricity the 
system is chaotic and the more the bubble is shifted from the center the 
larger part of phase-space is occupied by chaotic trajectories. Nevertheless, 
the shell effects are still strong. One can show that the shortest periodic 
orbit determines the gross structure of the shell energy when the bubble is 
close to the surface [9, 10]. 

The problem of two or more objects immersed in an infinite Fermi sys- 
tem has a similar character. In an infinite system the presence of impurities 
results in an appearance of resonances, which contribute to the shell cor- 
rection energy. In order to better appreciate the nature of the problem, let 
us consider the situation when two identical spherical bubbles have been 
formed in an otherwise homogeneous Fermi system. We shall ignore here 
the possible Coulomb interaction, as its main contribution is to the smooth 
part of the total energy of the system. In the semiclassical approach, which 
is justified for the "sizeable" bubbles (i.e. when the Fermi wavelength is 
small comparing to the size of the bubble), the shell correction energy is de- 
termined by the periodic orbits in the system. In the case of two spherical 
bubbles there exist only one such trajectory (with repetitions) which gives 
rise to the interaction energy between bubbles. It is the hyperbolic orbit 
lying in the line connecting the bubbles centers and characterized by the 

i + i + yji (l + 2 ) 



Lyapunov exponent: A = 2 In 



where R is the radius 



of the bubble and d is the distance between their centers. The interaction 
energy between the two bubbles due to the existence of this periodic orbit 
reads: 

Ti 2 kp 1 ^ [2nkF(lcos(2nkpd) — sin(2n/cp(f)] 
shell = — 47r{kFd)2 2_, n3 s inh 2 (nA/2) ' ( j 

where kF denotes the Fermi momentum and m the mass. Similar arguments 
can be presented for other shapes, e.g. cylinders or plates. One finds that 
at large separations the interaction energy oscillates as well but it decays as 
oc 1/d 5 / 2 in case of cylinders and as oc 1/d 2 in the case of plates. 

Hence it is clearly seen that there are some preferable arrangements 
in the systems of two impurities. The interaction between them depends 
mostly on their shapes and the geometry of the mutual arrangement. Our 
results show that shell effects associated with the appearance of inhomo- 
geneities in the neutron matter at subnuclear densities may strongly influ- 
ence the phase transition pattern in the neutron star crust since the shell 
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effects for various "nuclear lattices" are of the same order as the differences 
of the total energy between phasesfll]. 

We suspect that there are a lot of other effects, which might be rele- 
vant. We did not consider periodic orbits bouncing between three or more 
objects. An orbit bouncing between two bodies leads to a pairwise inter- 
action. Orbits bouncing between three or more bodies would lead to many 
body interactions, which however turn out to be of lesser importance [12]. 
We have also considered only perfectly smooth objects. If one allows for 
some degree of corrugation of these surfaces, many more periodic orbits are 
likely to appear and that would lead to even more complicated interactions 
and more complicated interference patterns. If temperature has a simple to 
predict qualitative effect on the shell correction energy, the role of pairing 
effects is not that transparent and deserves further investigation. 
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